J 


June  20, 1991  Riciiard  T.  Lih<-y,  Jr. 

lidmird  E.  Jr  I'ni'i  srur  nf  Ei:!^:>:r,  hul; 


Dr.  Edwin  P.  Rood 
Office  of  Naval  Research 
Scientific  Office  Code-1 132F 
800  North  Quincy  Street 
Arlington,  VA  22217-5000 

Dear  Dr.  Rood: 


The  purpose  of  this  letter  is  to  transmit  the  Second  Quarterly  Progress  Report  for  ONR  grant  N00014- 
91-J-1271  ("An  Experimental  Study  of  Plunging  Liquid  Jet  Induced  Air  Carryunder  and  Dispersion", 
Lahey  &  Drew  -  CoPI). 

The  second  quarter  of  this  research  project  has  been  devoted  to  completion  of  both  the  axisymmetric  and 
the  planar  nozzle  designs,  to  qualifjnng  the  DANTEC  3-D  LDA/PDA  system  for  underwater  opjeration, 
and  to  a  preliminary  analysis  of  two-phase  jet  spreading. 


EXPERIMENT 

As  discussed  in  th<»  last  quarterly  report,  an  axisymmetric  nozzle  will  be  used  to  benchmark  our 
experimental  results  against  prior  data.  Subsequently,  a  planar  nozzle  will  be  used  to  obtain  data 
which  should  be  more  representative  of  the  air  carryunder  process  induced  by  a  ship's  bow  wave. 

During  this  report  p>eriod  the  preliminary  nozzle  designs,  reported  on  last  quarter,  were  reviewed  by 
wind  tunnel  design  experts  at  RPI.  In  particular.  Professors  Nagamatsu  and  Hirsa  reviewed  our  nozzle 
design  and  based  on  their  comments,  modifications  were  made.  The  final  nuzzle  desigris  are  shown  in 
Figures  1-2  (axisymmetric  nozzle)  and  Figures  3-4  (planar  nozzle).  These  designs  have  now  been  turned 
into  the  RPI  machine  shop  and  will  be  fabricated  during  the  third  quarter  of  this  research  project. 

A  great  deal  of  exjjerimental  qualification  work  was  also  done  during  the  second  quarter.  In  particular, 
the  DANTEC  3-D  LDA/PDA  fiber  flow  system  was  designed  and  tested  underwater.  Unfortunately,  it 
was  found  that  the  sealing  of  the  lens  was  not  adequate  (ie,  leakage  occurred).  All  sets  of  lenses  (301 
mm,  600  nrvm  and  1200  mm)  have  been  returned  to  DANTEC  and  they  are  currently  being  redesigned  and 
manufactured  (in  Denmark)  to  obtain  a  leak-tight  seal.  Delivery  of  these  new  lenses  is  expect^  during 
the  third  quarter. 

The  Benjamin  System  3-D  traversing  mechanism  was  delivered  during  the  report  period  and  was 
installed  and  operated.  This  traversing  mechanism  has  the  capability  of  giving  one  micron  positioning 
accuracy  in  each  orthogonal  direction  (ie;  x,  y  and  z).  It  is  not  clear  at  this  point  whether  our 
measurements  must  be  made  in  the  forward  or  backward  scattering  direction.  If  forward  scattering  is 
required  then  a  special  mirror  system  will  be  needed  to  allow  for  accurate  PDA  pickup.  Such 
determinations  can  only  be  made  when  the  entire  system  is  in  place  and  operated.  In  accordance  with 
our  original  Work  Scope,  this  is  scheduled  for  next  year. 
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ANALYSIS 

Dxiring  the  Second  Quarter  of  this  research  project  preliminary  analysis  was  performed  on  two-phase 
jet  spreading  phenomena. 

Appropriate  equations  of  motion  have  been  derived  by  assembling  the  relevant  equations  from  Arnold's 
(1988)  PhD.  Thesis  (G.  Arnold,  "Entropy  and  Objectivity  as  Constraints  Upon  Constitutive  Equations  for 
Two-Fluid  Modeling  of  Multiphase  Hows"),  anid  adding  the  appropriate  term  representing  the  action 
of  surface  tension,  which  was  not,Jncluded,  in  Arnold’s  thesis.  These  equations  were  non- 
dimensionalized  for  a  long,  narrow  jet.  The  resulting  equations  are  parabolic,  so  that  numerical  solution 
should  be  straightforward.  - , 

In  addition,  we  have  approximated^ the  equations  of  motion  by  assuming  that  the  gas  volume  fraction  is 
small.  This  allows  us,To  use  the  known  solution  for  the  velocity  field  in  a  planar  jet  for  the  liquid.  If  we 
further  approximate  the  momentum  equation  for  the  bubbles  by  neglecting  ail  ertects  except  diag, 
virtual  mass,  and  liquid  pressure  gradient,  we  are  able  tp  solve  the  resulting  equations  for  bubble 
trajectories  in  the  known  liquid  velocity  field.  It  was  assunrted  that  the  bubbles  do  not  affect  the  liquid 
velocity  field  and  that  the  gas  density  is  negligible.  The  resulting  equation  in  Lagrangian  form  is 


=  ^g 


d  ci  1 

PgVb^Vg=p/Vb^v/(rg(t))  +  CvinP^Vb^(vi-Vg) +  2?/ AbCo  h-%  (Y/-Yg) 


+  ipt  -  pg)Vb  g 


(1) 


The  ToUmein  solution  for  a  planar  jet  was  used  for  the  liquid  velocity. 

The  system  (1)  was  solved  numerically  using  a  4*^  order  Runge-Kutta  subroutine.  The  computed 
trajectories  are  shown  in  Figure  5  for  different  initial  conditions.  Vg(t=0)  =  0  and  p  €  l-l,  1]  is  a 
parameter  that  specifies  rg(t=0).  See  Figure  6. 

One  result  of  this  calculation  is  the  point  s  where  the  trajectory  comes  to  the  surface.  Figure  7  depicts  s 
as  a  function  of  the  parameter  p.  It  can  be  easily  seen  that  for  p  near  -1  there  are  trajectories  that  cross 
each  other. 

If  P(p)  is  the  distribution  of  bubbles  at  the  inlet  and  S{s)  is  the  distribution  at  the  free  surface,  then 


S(s)  = 


dpi 

ds 


P[pi(s)]  + 


dp2 

ds 


P[p2(s)] 


where  p^(s)  and  P2(s)  are  nwnotone. 
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Improvements  to  this  model  are  contemplated,  including  evaluation  of  the  effects  of  the  neglected 
terms,  and  using  this  model  with  a  realistic  bubble  density  distribution  at  the  beginning  of  the  jet,  in 
order  to  obtain  tf»e  bubble  density  field. 


This  last  concept  points  to  a  research  direction  that  we  shall  pursue;  namely,  the  developnwnt  of  the 
aur  being  dragged  along  the  liquid  jet  into  bubbles  under  the  surface.  We  have  been  developing  a  model 
for  bubble  coalescence  and  breakup  for  predicting  flow  regime  transition  for  bubbly  flow.  We  propose  to 
use  this  model  to  see  if  it  is  capxable  of  predicting  bubble  formation  for  the  plunging  jet. 


If  you  have  any  questions  concerning  the  material  presented  in  this  letter  please  do  not  hesitate  to 
contact  us  [Lahey;  518-276-8579;  Drew:  518-276-6903]. 


Sincerely  yours,  , 

kliiUi  0\  t 

Di.  D.A.  Drew 

Professor  of  Mathematical  Sciences 


DAD/ev 


cc 


Administrative  Grants  Officer  -  ONR 
Director  -  Naval  Research  Laboratory 
Defense  Technical  Information  Center 
R.T.  Uhey,  Jr.  (RPI) 

F.  Bonetto  (RPI) 
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FIGURE-2  INLET  (CALMING)  SECTION  FOR  AXISYMMETRIC  NOZZLE 
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FIGURE-4  INLET  (CALMING)  SECTION  FOR  PLANAR  NOZZLE 
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